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ABSTRACT 

We use cosniological SPH simulations to study the effects of mergers in tlie star forma- 
tion history of galactic objects in hierarchical clustering scenarios. We find that during 
some merger events, gaseous discs can experience two starbursts: the first one during 
the orbital decay phase, due to gas inflows driven as the satellite approaches, and the 
second one, when the two baryonic clumps collide. A trend for these first induced star- 
bursts to be more efficient at transforming the gas into stars is also found. We detect 
that systems which do not experience early gas inflows have well-formed stellar bulges 
and more concentrated potential wells, which seem to be responsible for preventing 
further gas inward transport triggered by tidal forces. The potential wells concentrate 
due to the accumulation of baryons in the central regions and of dark matter as the 
result of the pulling in by baryons. The coupled evolution of the dark matter and 
baryons would lead to an evolutionary sequence during which systems with shallower 
total potential wells suffer early gas inflows during the orbital decay phase that help 
to feed their central mass concentration, pulling in dark matter and contributing to 
build up more stable systems. Within this scenario, starbursts triggered by early gas 
inflows are more likely to occur at early stages of evolution of the systems and to 
be an important contributor to the formation of stellar bulges. Our results constitute 
the first proof that bulges can form as the product of collapse, collisions and secular 
evolution in a cosmological framework, and they are consistent with a rejuvenation of 
the stellar population in bulges at intermediate z with, at least, 50% of the stars (in 
SCDM) being formed at high z. 
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1 INTRODUCTION 

The physical mechanisms that trigger and regulate star for- 
mation in galaxies are thought to be of different nature and 
act at different scales (e.g., Kennicutt 1998). Among them, 
mergers and interactions are known to play a crucial role as 
it has been demonstrated from observations of high emis- 
sion levels of Ha (e.g., Kennicutt et al. 1987), far infrared 
(e.g., Solomon & Sage 1988; Sanders & Mirabel 1996) and 
radio continuum (e.g., Condon et al. 1982) in interacting 
systems. Recent observations of nearby galaxy pairs by Bar- 
ton, Geller & Kenyon (2000; see also Donzelli & Pastoriza 
1997) found that interactions during close encounters can 
be statistically correlated with enhancements of the star for- 
mation activity lasting more than 10* yr. It has also been 
observed that moderately luminous starbursts in the nearby 
Universe may occur in disc galaxies that show only mild 



external perturbations originated in moderate interactions 
(Gallagher et al. 2001). Regarding our Galaxy, a new study 
of the chemical composition of the solar neighbourhood by 
Rocha-Pinto, Maciel & Flynn (2000) supports the possible 
existence of three main starburst episodes, lasting each one 
few Gyr, which can be associated with nearby passage of 
the Milky Way galaxy satellites. Remarkably, two of such 
starbursts show evidences of being associated to only one of 
those close encounters. Recent observations of high redshift 
(z) objects show an increase of morphologically disturbed 
and interacting objects with z, accompanied by an increase 
of the star formation activity (e.g., Bouwer, Broadhurst & 
Silk 1998; Driver et al. 1998; Brinchmann et al. 1998; Saw- 
icki & Yee 1998; Le Fevre et al. 2000). Ahhough it is not 
yet clear at what extent these data are fully consistent with 
predictions of models based on hierarchical clustering, all 
of them agree in supporting a scenario in which interac- 
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tions/mergers and star formation increase with z (Flores et 
al. 1999; Steidel et al. 1999). 

Understanding the nature of galaxies as a function of 
z and the interface between galaxy formation and structure 
formation, requires to be able to describe how star formation 
proceeds in galactic objects at different stages of evolution, 
and how this activity is affected by the formation of the 
structure at larger scales. This is a complex task taking into 
account that as a galaxy evolves, it may react to the same 
physical mechanism (i.e., such as mergers/interactions) in a 
different way, depending on its astrophysical and dynamical 
properties (i.e., formation of a bulge, baryonic/dark matter 
dominated potential well at the centre, availability of gas, 
etc), and, consequently, on its previous evolutionary history. 

In a hierarchical clustering scenario, galactic objects 
form by the aggregation of substructure. Assuming this 
model as the more feasible one involves the understanding 
of the role played by mergers and interactions as an object 
evolves. In particular, pre-prepared simulations of merging 
galaxies show that, during violent events, strong gas inflows 
can be triggered, feeding star formation activity at the cen- 
tre (e.g., Noguchi & Ishibashi 1986; Noguchi 1991; Mihos, 
Bothum & Richstone 1993; Mihos 1992; Barnes & Hernquist 
1991, 1996). Furthermore, Mihos & Hernquist (1994, 1996 
hereafter MH96) demonstrated that during merger events, a 
bulge-less disc galaxy is more susceptible to growing bar in- 
stabilities that might trigger gas inflows towards its centre, 
even before the actual collision of the baryonic clumps oc- 
curs, than a bulge-disc system. In the first case, a double star 
formation burst was identified, with a first component trig- 
gered by the strong tidal fields developed during the orbital 
decay phase (hereafter ODP), and a second one triggered 
when the baryonic cores collided. The relative strength of 
the double bursts was shown to depend on the compactness 
of the bulge and on the availability of gas at the time of the 
merger. Conversely, a disc galaxy with a stellar bulge was 
found to be more stable, and only when the actual collision 
of the baryonic cores occurred, a starburst was produced. 
These authors directly link the response of the systems dur- 
ing a merger to the characteristics of their potential wells. 

However, iu these pre-prepared mergers, the great va- 
riety of phenomena involved in the complex processes of 
galaxy-like objects (GLOs) building-up from the field of pri- 
mordial fluctuations at high z up to z = 0, cannot be fol- 
lowed. Such phenomena include collapse events, mergers and 
interactions with close neighbours, accretion of smaller ob- 
jects, gas infall at scales of some hundreds of kpcs (from 
the environment to the halo or from the halo to the disc), 
gas inflow at scales of tens of kpcs (within discs and from 
discs to the bulge-like central concentrations), etc. All of 
them play a role in the assembly of galaxy-like objects, and, 
moreover, they interact among themselves in a non-trivial 
way. Pre-prepared mergers are very detailed experiments, 
but which only describe a particular phase of galaxy forma- 
tion from initial conditions that are generally set by hand. 
Moreover, mergers in these experiments are taken in isola- 
tion, and not in connection with the other processes involved 
in galaxy formation. On the contrary, with a self-consistent 
approach (i.e., following the assembly of galaxy-like objects 
from the fleld of primordial fluctuations at high z up to 
« = 0, in the framework of a cosmological model), mergers 
can be studied in connection with other mechanisms, that 



play a fundamental role in the merger process itself, such as 
gas infall at some hundreds of kpc scales and gas inflows in 
discs. Also, one can analyze how the initial conditions that 
are usually considered in pre-prepared simulations (at scales 
of about their virial radii, i.e., some hundreds of kpcs), ap- 
pear, along the evolution of the objects, in connection with 
the evolution of their environment at larger scales. At these 
larger scales, evolution is determined by the global evolu- 
tion of the cosmological model. Of particular interest in this 
regard is to learn how bulges are assembled, and how the 
growth of a given bulge at a given z might prevent it from 
growing at lower zs. 

By using cosmological hydrodynamical simulations, 
Dommguez-Tenreiro, Tissera & Saiz (1998, hereafter 
DTTS98) and Tissera, Saiz & Dommguez-Tenreiro (in 
preparation) analyzed how gas inflows can be generated dur- 
ing the evolution of a typical galaxy-like object in a hierar- 
chical clustering model, as it forms by the merger of sub- 
structure. These authors showed that the gas with non-zero 
angular momentum always tends to form a disc, but the 
survival probability of these discs depend on their stability 
properties. Their findings confirmed that gas inflows can be 
easily triggered during violent events if the axi-symmetrical 
character of the potential well is not assured by a central 
mass concentration, as for example, a dark matter halo dy- 
namically dominating at the centre, or by the presence of 
a compact stellar bulge (e.g.,Toomre &i Toomre 1972; Mar- 
tinet 1995 and references therein). These gas inflows can be 
mild or catastrophic, depending on the internal properties of 
the disc-like systems. Catastrophic inflows can destroy the 
discs and are easily induced in the absence of a compact 
central mass concentration (DTTS98; Steinmetz & Navarro 
1999). However, if stellar bulges are allowed to form, discs 
with observational counterparts can be formed in hierarchi- 
cal clustering scenarios (Saiz et al. 2001, hereafter SalOl). 

Tissera (2000a, hereafter TisOO) studied the star for- 
mation history of galactic-like objects in relation to their 
merger histories within a cosmological context. This author 
found that minor and major mergers may trigger staxbursts 
as powerful, regardless of their relative masses, if there is 
enough available gas in the system. Only the starbursts ef- 
flciencies, defined as the fraction of the available gaseous 
mass in the system that is actually transformed into stars, 
were found to depend on the relative masses of the colliding 
objects. The stellar masses formed and starbursts (SBs) du- 
rations measured in the simulations are in agreement with 
observational results (e.g.. Barton et al. 2000; Glazebrook et 
al. 1998). Even more, TisOO estimated a change of up to 0.6 
magnitudes in the colours of the simulated galaxies due to 
mergers, consistent with observational results reported by 
Le Fcvre et al. (2000). It was also pointed out by the author 
that during certain merger events, double SBs were present. 
If these phenomena were a common physical process, then, 
understanding how they could be induced in a hierarchical 
clustering scenario is very relevant to the study of galaxy for- 
mation. The physical interpretation of several observations 
of galaxy properties, such as luminosity function or spatial 
distributions, may be affected by the possibility of induc- 
ing nmltiple starbursts during a single merger event (Tis- 
sera 2000b). Moreover, the formation of the stellar bulges of 
spiral galaxies can be directly affected by the possibility of 
combining two feasible formation mechanisms, such as sec- 
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ular evolution and collisions, during a single merger event. 
Observationally, there are evidences that suggest that these 
different processes, together with collapse at high z, could 
work with different relative importance, depending on the 
Hubble type, to manufacture stellar bulges (Ellis, Abraham 
& Dickinson 2001). 

In this paper, we focus on the study of double star- 
bursts within the context of galaxy formation in a cosmo- 
logical framework. Contrary to pre-preparcd mergers, we 
use fully solf-consistont cosmological hydrodyiiaiiiical sim- 
ulations, where the distribution of merger parameters, such 
as orbital characteristics, the orbital energy and angular mo- 
mentum, the masses of the virial haloes and baryonic clumps 
involved in the merger event, and the spin, internal struc- 
ture, and relative orientation of the baryonic clumps that 
are about to merge, among others, arise naturally at each 
epoch as a consequence of the initial spectrum of the den- 
sity fluctuation field, its normalization, and the cosmological 
models and its parameters. The effects of mergers and in- 
teractions on star formation can be studied in a consistent 
scenario, although at the expense of losing numerical res- 
olution. In Section 2 we present the main aspects of the 
numerical models and in Section 3, 4 and 5, the analysis of 
starbursts. Section 6 summarizes the results. 



2 NUMERICAL MODELS 

The simulations analyzed take into account the gravitational 
and hydrodynamical evolution of the matter including an 
algorithm to transform the cold and dense gas into stars. 
We use the Smooth Particle Hydrodynamical (SPH) code 
developed by Tisscra, Lambas & Abadi (1997) coupled to 
the AP3M (Thomas & Couchman 1992). SPH based codes 
have proved to be powerful tools to study galaxy formation 
(e.g., Navarro & White 1994; Steinmetz & Navarro 1999; 
DTTS98; TisOO; SalOl). 

We run four simulations, S.l, S.2, S.3 and S.4, consistent 
with Cold Dark Matter universes with Q = 1, A = 0. The 
normalization of the power spectrum was taken to be as = 
0.40 for S.l and S.2 and ag = 0.67 for S.3 and S.4 (see Ta- 
ble 1). The simulated boxes have 5h~^ Mpc length {h = 0.5) 
with N = 64* total particles. Baryonic particles represent 
10 % of the total mass. Note that dark matter and bary- 
onic particles have the same mass, Mpart = 2.6 x 10* M0. 
The gravitational softening used is 3 kpc, and the smaller 
smoothing length allowed is 1.5 kpc. Simulations S.l and 
S.2 share the same initial conditions (A) while S.3 and S.4 
are different realizations of the power spectrum (B, C). Run 
S.l corresponds to experiments S.l in DTTS98, SalOl and 
TisOO. Simulation S.2 has been run with a higher star for- 
mation efficiency in order to assess the dependence of the 
results on this parameter. 

The star formation (SF) algorithm transforms dense 
and cold gas in a convergent flow into stars according to 
the Schmidt law (Navarro & White 1994; TisOO). We as- 
sume a certain star formation efficiency, c (see Table 1), 
and a characteristic time-scale for the star formation pro- 
cess. The SF parameters used in S.l, S.3 and S.4 allow the 
transformation of gas into stars only in the very dense re- 
gions. In practice, this condition implies that gas particles 
have to fall in all the way down to the very central regions 



of GLOs before they arc converted into stars (see Fig. 2). 
This assures the formation of a compact central stellar mass 
concentration, that is, a bulge, and that the transformation 
of gas into stars on the discs is suppressed. Conversely, the 
SF process in S.2 is more efficient, which implies that gas 
particles are transformed into stars sooner than those in S.l. 
This difference translates in different physical properties for 
the final objects, as Tissera & Dommguez-Tenreiro (1998) 
have shown. These authors study the mass distribution of 
galactic objects in S.l and S.2 finding that, because of the 
higher star formation efficiency, those in S.2 tend to be dom- 
inated by stellar spheroidal components with potential wells 
not as deep as those of their counterparts in S.l. The results 
of S.2 have been included in this work to show that the con- 
clusions we reach here are not mainly determined by the 
star formation parameters, implying that the mechanism at 
work is a general physical process. 

We want to point out that, because the SF process is 
based on the Schmidt law, any sudden increase of the gas 
density directly implies a corresponding increase in the star 
formation activity. According to TisOO (see also MH96) , it is 
possible to have a large percentage of gas that satisfies the 
SF criterion, but that is transformed into stars at a quiescent 
rate. And, on the other hand, a smaller fraction of gas can 
be suddenly compressed, triggering a starburst. 

Supernova feedback effects have not been included in 
these simulations. Energy feedback processes are believed 
to play a key role in helping to set a self-regulated star for- 
mation regimen (Silk 2001). However, these processes are 
not well understood and its modelling in hydrodynamical 
simulations is still quite controversial (Katz 1992; Navarro 
& White 1994; Metzler & Evrard 1994; Yepes et al. 1997; 
Moscoru et al. 2001). And so, as a first step, it is wise to 
analyze stax formation separately from supernova effects. 
Furthermore, as pointed out by DTTS98 and SalOl, the re- 
strictive SF model used in this work could have allowed us, 
in a way, to mimic the effects of supernova energy feedback. 
In fact, the SF algorithm wo have used allows the formation 
of stellar bulges without quickly depleting the gas reservoirs 
of the objects. SalOl have shown that the formation of com- 
pact stellar bulges provides stability to the disc-like objects, 
so that the disc component is able to conserve an important 
fraction of its angular momentum during violent events, pro- 
ducing systems that resemble current normal spirals. Partic- 
ularly, these authors have analyzed the dynamical and struc- 
tural properties of GLOs in S.l at z — 0, finding that they 
have observational counterparts. This is also valid for those 
GLOs in S.3 and S.4 that resemble disc-like objects. Hence, 
we are working with a set of objects that can be fairly com- 
pared to spiral galaxies at z = 0. 



3 STAR FORMATION HISTORIES AND 
MERGERS 

Galaxy-like objects are identified at their virial radius, r2oo 
(i.e., the radius where 5p/p « 200; White & Frenk 1991) 
at z = 0. We reject those GLOs with a massive comparuon 
within two virial radii at z = to avoid complications due to 
tidal fields originated by the presence of companions and/or 
the underlying overdensity, focusing on the effects produced 
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by the assembly of each individual object through hierar- 
chical growth. We are only going to study those GLOs with 
more than 250 baryonic particles inside the main object, 
and, generally, from z = 1 to diminish as much as possi- 
ble numerical resolution effects, which could be present at 
higher z. No other constraints have been imposed to select 
the analyzed GLOs. 

A typical GLO is composed by a dark matter halo, a 
main baryonic clump and a series of small baryonic satellites. 
Table 2 gives the total number of dark matter (A^dark) and 
baryonic (Abar) particles within the virial radius for each 
GLO at 2 = 0. Note that dark matter haloes are very well 
resolved by several thousand particles. According to results 
from Steinmetz & White (1997), if the potential well of a 
galactic halo is well-reproduced, the dissipative component 
is forced to set onto the correct density profile, even if it is 
resolved by few hundred particles. If the gas density profile is 
adequately described, then the SF process within the GLO 
can be rehably followed (see also DTTS98). 

The merger trees of the GLOs are reconstructed by 
tracking back all particles that belong to a GLO at z = 0. 
The progenitor object is chosen as the more massive clump 
within this merger tree. A merger event will be defined as the 
whole process from the instant when the two main baryonic 
objects are first identified to share the same dark matter halo 
(this instant is the beginning of the ODP, whose redshift will 
hereafter be termed Zo), until the instant they actually col- 
lide at a redshift (once the ODP is over). Note that, in 
this paper, a collision means that the baryonic systems can- 
not be longer separated into two distinctive objects, but it 
does not imply any particular characteristics for the rem- 
nant. During the ODP, the minor colliding baryonic clump 
will be referred to as the satellite. 

The corresponding star formation rate (SFR) histories 
of the progenitors can be constructed as the ratio between 
the stellar mass formed during a time-step of integration and 
its duration (At = 1.3 x 10^ yrs for S.l and S.2 and At = 
1.2 X 10^ yrs for S.3 and S.4). These SFR histories can be 
described as a result of two contributions, a quiescent SF and 
a series of starbursts (see TisOO for details). The quiescent 
SF gives rise to a constant SFR or threshold; starbursts are 
defined to be those points of the SFR history of a given 
object above its threshold (see Fig. 1). 

Mergers of substructure are responsible of the triggering 
of some of these SBs. ^From all SFR peaks identified, in this 
paper we are only concerned with those that occur within 
merger events. We have studied 29 merger events related to 
43 starbursts of which 14 are triggered during orbital decay 
phases. Once the starbursts are isolated, it is direct to define 
their maximum SFR, (Tstar, the total stellar mass formed, 
Mburst, and their duration, Tburst- These quantities deter- 
mine the characteristics of each starburst, together with the 
gas richness of the system, Mstar/Mbar, and the virial mass 
ratio of the colliding objects, Msat/Mpro, at the time of the 
merger. Mstar and Mhai are the stellar and baryonic (i.e., gas 
plus stars) content of the merging systems within r2oo; Msat 
and Mpro are the virial masses of the satellite and progeni- 
tor, respectively. Table 2 gives their values for the mergers 
analyzed in this work. 



4 DOUBLE STARBURSTS 

During some merger events, two starbursts can be identified 
in the SFR histories of the GLOs. We classify as a double 
SB those two well-defined SFR peaks that could be directly 
related to the same merger process. In other words, those 
two bursts that appear in the time interval from Zo (i.e., 
the instant when the satellite enters the virial radius of the 
progenitor) to the actual collision of the baryonic cores at 
Zc ■ Only those bursts that strictly satisfy this condition were 
taken as probably induced by the same merger. Whenever 
a doubt exists, we choose to eliminate the merger from our 
statistics so that it does not contribute as either, a single 
burst or a double one. This situation, actually, occurs once. 
As an example, we show in Fig. 1 the SFR for typical GLOs 
that experience (a) a double SB (GLO 3 in S.l, Fig. la) 
and (b) a single one (GLO 3 in S.3, second merger. Fig. lb) 
as a function of lookback time {t{z) = 1 — (1 + z)~^^'^). 
We also plot in each of these Figures the distances from the 
progenitor to the satellite in the time interval from Zo to 
ZcQ. In Fig. la, we see that the first burst occurs during the 
orbital decay phase, while the second one is triggered during 
the collision of the baryonic cores. In Fig. lb, only at the 
impact of the two baryonic clumps an enhancement of the 
SF activity is detected. Note that when the satellite enters 
the virial radius of the progenitor, in both cases (a, b), the 
progenitors are equally gas-rich as can be appreciated from 
Table 2. 

For each double SB, we measured the rbmst, Afburst, 
(Jstar and the time interval Atburst between the primary 
(""star) and the secondary (o"star) burst maxima, (zburst.i > 
Zburst,2; Atburst was mcasured from the last SFR value of the 
primary burst to the first point of the secondary one, above 
the thresholds defining the bursts). The temporal separa- 
tions between the bursts, Atburst ~ 5 x 10* — 5 x 10^ yrs, 
are consistent with the collapse times of satellites estimated 
from dynamical friction considerations (Navarro, Frenk & 
White 1995; Hernquist 1989). Table 2 summarizes the re- 
sults (secondary components of double bursts are indicated 
by the letter D). Last column on Table 2 gives the redshifts 
[z) at which each starburst starts. 

It is clear that if there is enough available gas in the sys- 
tem, a starburst will be triggered when the progenitor and 
the satellite collide at z^. Regarding double SB events, our 
concern now is to prove that the SBs detected in the progen- 
itors during some ODPs are produced by gas inflows driven 
as the satellite gets closer, as predicted by pre-prepared sim- 
ulations (e.g., Hernquist 1989; MH96). Also, we want now 
to assess the differences, prior to the merger event, of the 
internal structure of GLOs involved in double and single 
SBs, as well as the evolution of these differences, in relation 
with their SFR histories, during and immediately after the 
merger event. To this end, we have studied and compared 
the properties of our GLO sample at the merger events re- 
ported in Table 2. For each of these mergers, the variations 
of the baryonic mass and angular momentum distribution 
at the central regions of their respective progenitors have 
been followed and quantified, in the time interval between 

* It is worth mentioning that in some merger events, during the 
orbiting of the satellite around the progenitor, more than one SBs 
can be detected. These cases are not studied in this paper. 
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the beginning of their ODPs at Zo and a reference time at 
Zi^, using different estimators. 

As a first step, the progenitor particle number (actu- 
ally, the baryon mass) inside a given radius, A''(< r), has 
been measured at different redshifts in the [2:o,2:f] time in- 
terval. The general result is that gas inflows are detected 
during this time interval at the central regions of progeni- 
tors that experience double starbursts. In contrast with this 
result, no gas inflows are detected at any stage at the cen- 
tral regions of progenitors involved in merger events with 
single starbursts. To illustrate these results, we have chosen 
as typical examples the mergers occurring in GLO 2 from 
S.4 (a double SB event a.t z ^ 1.20) and in GLO 5 from S.l 
(the latest one at z « 0.19, a single SB event; see Table 2), 
but the behaviour is similar for the other merger events. In 
Fig. 2a we plot A^(< r) for the primary component of the 
double SB in GLO 2 from Zo = 1.81 to zi = 1.28. Note how 
gas particles are displaced from the outer to the inner ra- 
dius as the time goes on. Part of the particles that get into 
the 10 kpc < r < 30 kpc shell are transported all the way 
down to the very centre (r < 4 kpc). This increase in the 
gas density immediately translates into an increase of the 
star formation activity, as previously discussed. The inward 
displacement of these particles is due to a loss of their spin 
angular momentum content as the satellite gets closer (see 
DTTS98). GLO 5 has assembled most of its mass at higher 
z, so that when the last merger occurs, its gaseous disc and 
stellar bulge are clearly at place. Fig. 2b shows A'^(< r) be- 
tween Zo and Zf for the last merger. As can be seen, the 
changes in the baryonic mass distribution are much less im- 
portant in this case than for GLO 2. There is almost no 
increase in the baryonic mass within 15 kpc or the inner 4 
kpc, implying that no significant gas infiows have been trig- 
gered in the central region during the orbital decay phase of 
this satellite. 

In the other six plots of Fig. 2, we show the angular mo- 
mentum component per unit mass, jz,i, along the direction 
of the total angular momentum of the main baryonic sys- 
tems, for each gas and star particle in GLO 2 (left panels) 
and GLO 5 (right panels), versus particle radial distance, 
ri, at three stages of evolution within the interval of interest 
[zo, Z{]. These plots give a clear picture of the internal struc- 
ture of both GLOs during such interval. These Figures show 
that when GLOs have well-formed stellar bulges, the disc 
components are not significantly affected during the ODP. 
Conversely, those GLOs which experience early gas inflows 
do so at the expense of the gas component on the discs. 
These gas inflows feed the stellar bulges. Fig. 2 also sug- 
gests that in the latter cases, GLOs are in early stages of 
evolution. 

To go further into the quantification of the differences 
between the properties of the GLOs at a double or single 
SB event, we use results from bulge-disc decomposition of 
their projected mass surface density (Scannapieco et al., in 



t For mergers producing double starbursts, Z{ has been taken 
to be the moment when the ODP ends, 2o,fi because [20,^0 f] is 
the time interval within which the primary component of double 
starbursts is detected. For mergers that do not produce double 
starbursts, we have taken Zf = Zc, that is, the time when the 
satellite and the progenitor actually collide. 



preparation). Following SalOl, the projected mass surface 
density has been decomposed by using a Sersic law for the 
bulge component and an exponential law for the disc compo- 
nent. SalOl have shown that, at 2: = 0, this procedure yields 
scalelengths and dynamical parameters which are in very 
good agreement with recent observations of spiral galaxies 
(e.g., Courteau 1996). Here, the aim of this decomposition 
is to provide structural parameters which allow a convenient 
quantification of the mass distributions of progenitors at Zo 
and of their changes between Zo and zf. The analysis has 
been performed for GLOs in S.l and S.4, for which a suit- 
able number of outputs to follow their evolution in time are 
available. The general result is that the mean values of the 
masses of the bulges at Zo, as well as their variations between 
Zo and 20,5, show clear differences for GLOs involved in dou- 
ble or single starburst events. At Zo, the former show smaller 
bulges than the later on average (i.e., « 50% smaller), but, 
after the starbursts that take place within the [2o,Zo,f] in- 
terval, their bulges have grown by 47 ± 19% on average. 
Meanwhile, those GLOs with no important SF during an 
equivalent time-interval, show a change in the bulge mass 
of 5 ± 2%0 These results confirm that systems experiencing 
early gas inflows have smaller bulges which are fed by gas 
infiows during the ODPs. 



5 STARBURST PROPERTIES 

Taking into account previous works on disc formation, there 
are evidences that support the idea that the stability of the 
disc structures plays a relevant role in the triggering of gas 
infiows when the systems are subject to tidal fields. This 
stability could be assured by a compact stellar bulge when 
the dark matter (DM) halo is not dynamically dominating at 
the centre (MH96; Barnes & Hernquist 1996). Then, a first 
step towards the understanding of the physical origin of the 
induced SBs in our simulations, is to analyze any possible 
relation between the formation of the stellar bulge in the 
progenitor and the burst properties. 

To this end, our first task is to define the epoch of bulge 
formation. This is a controversial point since bulges can con- 
tinuously grow by mergers, bar-driven infall or other mecha- 
nisms (Combes 2000). The simplest possibility to define the 
presence of a bulge is to adopt a threshold criterion relative 
to the bulge mass. However, we want to employ a criterion 
which allows the comparison among bulges of different virial 
masses. Consequently, due to the correlation among bulge 
masses and virial masses, a threshold criterion relative to 
normalized bulge masses should be used instead. Moreover, 
the normalizing factor has also to be sensitive to the history 
of evolution of GLOs (i.e., their respective merger trees and 
interaction events), because this history also affects the SF 
history of each GLO. We found that the total stellar mass 
in the GLO at z = is sensitive to both the virial mass 
of the object and its history of evolution. Then, we define 
the ratio between the stellar mass formed in the progenitor 
at a given 2, Afg^tar, and the total final stellar mass, A^fta^^", 

t The errors correspond to standard deviations calculated by 
using the re-sampling bootstrap technique with 500 random 
samples. 
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of the GLO, as an estimator of the stellar core (i.e., bulge) 
formation (see Fig. 2 for a general picture of baryon distri- 
butions). The monitoring of this ratio along the history of 
evolution of a GLO allows us to adopt a clean criterion for 
the formation of bulges. 

Now, if our hypothesis about the ability of well-formed 
bulges to provide stability to a bulge-disc system is correct, 
then there should be a correlation between the Mg^tar/A/star" 
ratios, measured at a redshift z = Zo, when the satellite en- 
ters the virial radius of the parent object, and the properties 
of the starbursts. For each merger in Table 2, we have cal- 
culated Ms^tar/A'/sta° and the relative strengths of the two 
bursts, o-^tiii / '^'stiii ■ The results arc plotted in Fig. 3. As can 
be seen from this Figure, there is a clear segregation which 
indicates that when a single burst ((Tstar/fJtar = 0) is associ- 
ated with a merger, GLOs tend to have larger M^^^^/M^^^l? 
values. Conversely, when a double burst is found, all of them 
have Mfftaj/Mg^tar' ^ 0.40. Our results show that as the 
stellar bulge forms, the relative strength relation reverses 
(""star/fJtar < 1) Until the first burst disappears when the 
stellar bulge is fully present (o-gtar/c'star = 0). The ratio 
Mrtar/Mfta° ~ 0.40 seems to indicate when the stellar bulges 
have grown large enough to prevent oven mild inflows. 

Hence, it seems plausible that, depending on the struc- 
tural properties of the GLOs, the merger with a satellite may 
induce two starbursts whose properties would also depend 
on the internal structure of the galactic objects. Our find- 
ings also agree with Hornquist (1989) results, which suggest 
that tidal triggering of starbursts may be a natural way of 
explaining the discrepancy found between the inferred star- 
burst time-scales (^ 10* yr) and the dynamical time-scales 
of the objects (^ lO^yr) . 

5.1 Dependence on the potential well 

In order to intend to relate the characteristics of the poten- 
tial well in the very central regions of GLOs with the pres- 
ence of a well- formed bulge as measured by Mgts,! / M^i^^ , we 
estimate the circular velocity at r = 2 kpc, w^, for the dark 
matter {k = DM) and gaseous (k = gas) components and 
the total mass (fe = tot) at Zo- This is a simple but straight- 
forward way of measuring the mass concentration. For this 
purpose, we calculate the circular velocity curves defined 
as V{r)^ = GM{r)/r, for the three mass components at 
primary components and single bursts. We use velocities in- 
stead of density profiles since the former are defined by the 
integrated mass which is less affected by numerical noise. 
Note also that for this calculation we are not including the 
gravitational softening since we are interested in measuring 
the actual mass distribution. However, the correct circular 
velocity curves should be estimated by taking into account 
such softening as shown by Tissera & Dommguez-Tenreiro 
(1998) and SalOl. 

Figs. 4 (a,b,c) show for the total mass, DM and 
gaseous components, respectively, against the ratio between 
the heights of the starbursts (Tstar/fstar (recall that for single 
bursts, cTstar/cstar = 0) for bursts in S.l and S.4. The segre- 
gation is quite clear, indicating that only when the systems 
have less concentrated potential wells measured by lower Vk 
values, gas inflows are induced, triggering starbursts dur- 
ing the orbital decay period. Note that all mass components 
show a similar trend. As it can be appreciated from Fig. 3, 



systems with o-gtar/fstar 7^ have smaller stellar bulges. Be- 
cause of the new gas infall, the stellar bulges grow as the 
gas is transformed into stars at the centre. When GLOs have 
well-formed bulges, they have already pulled in enough mass 
to assure the stability of the system which prevent further 
important inflows. 

5.2 Starburst efHciencies 

^From Fig. 3 we showed that double bursts occur only when 
Mit^r/M^{l° < 0.40. It was found by TisOO that the effi- 
ciency in the transformation of gas into stars during SBs 
correlates with Msat/Afpro in the sense that massive merg- 
ers (Msat/Afpro > 0.40) seem to be more efficient at convert- 
ing gas into stars. Here we explore the possibility that the 
SB efficiencies (defined as the ratio between the stellar mass 
formed during a burst and the total gas mass available in the 
system at the time of the merger, Mburst/Afgas) of the bursts 
triggered during the ODP may be different from those of SBs 
produced at the actual collision of the baryonic clumps. In 
double SBs we do not consider the secondary components, 
since they form from the left-over gas in the progenitor af- 
ter the primary one and the gas component associated to 
the satellite. They are also difficult to study since we would 
need to temporary resolve the merger with a higher number 
of outputs and higher numerical resolution. 

Fig. 5 shows Mburst/Mgas vs Mi^^JM^^Z"- Among GLOs 
in the same simulation, the larger values arc mainly recorded 
for the smaller M^t^^ / M^^^ . This can be clearly appreciated 
from Table 3 which shows the mean values of Mburst/.^gas 
for Miftar /-^^s^tar^ Smaller or greater than 0.40. Note that mean 
starburst efficiencies are higher for GLOs in S.2 simulation. 

Hence, we find a clear trend for SBs induced during 
the ODP to be more efficient at converting gas into stars 
than those produced by the actual collision of the baryonic 
cores. This result may have strong implications for galaxy 
formation. 

5.3 Dependence on redshift 

In Figs. 6a and 6b wc plot Msax/Mpro and M^^^^^ / M^uS , re- 
spectively, as a function of the redshift when the starbursts 
start (2:burBt on Table 2). We include single bursts and pri- 
mary components of double SBs (features with arrows). As 
can be seen form Fig. 6a, there is a clear trend to have major 
mergers (i.e., similar mass objects) at high z as expected in a 
hierarchical clustering model (e.g., Lacey & Cole 1993). We 
also see that double SBs tend to occur at higher z, and con- 
sequently, to be linked to more important merger events, 
with Msat/Mpro > 0.20 for aU double events. In Fig. 6b, 
A£ftar/A^s1;ar'' BiS & fuuctiou of Zburst shows how, as the pro- 
genitors build up their stellar bulges as they evolve, the 
probability of experiencing bursts induced via tidal forces 
diminishes. Note that, in these models, all double events 
occur at z > 0.4. 

^From the combination of the information coming out 
from Fig. 6, we note that 61% of bursts induced during the 
ODP are triggered by minor mergers (Msat/Afpro < 0.40). 
Even more, if we consider only SBs in S.3 and S.4, which 
constitute an homogeneous sample (i.e., different random 
phases of the power spectrum run with the same model pa- 
rameters), the trend is very similar: 66% of the SBs detected 
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during the ODP are related to minor mergers. This tendency 
suggests that the relative mass of the colliding objects is not 
the determining mechanism in the triggering of the early gas 
inflows. On the other hand, none of the systems with well- 
formed stellar bulges, M^^^^/M^^° > 0.40, have starbursts 
induced during the ODPs (Fig. 3). 

i,Prom these results, we conclude that hierarchical clus- 
tering is determining the global trend for massive mergers 
to occur at higher z, when the systems are more likely to be 
unstable, and, consequently, more susceptible to experienc- 
ing gas inflows and double SB events. We predict that this 
trend to have SBs induced via tidal forces at higher z, com- 
bined with their higher efficiency in transforming gas into 
stars, lead to a global star formation efficiency that would 
increase with z. Observational analysis of the star formation 
efficiency with z should test this prediction. 

5.4 Bulge formation 

Three different formation scenarios have been proposed 
for bulges of spiral galaxies: dissipative collapse (Eggen, 
Lynden-Bell & Sandage 1962), mergers of previously formed 
systems (e.g., Kauffmann, White & Guiderdoni 1993) and 
secular evolution (see Combes 2000 for a review). It has 
been also suggested that the three of them may have played 
a role, but with different relative efficiencies depending on 
the Hubble type (e.g.. Combes 2000; Elhs ot al. 2001). Our 
results support such hypotheses. Moreover, tlioy show that 
a given merger event could contribute to the formation of 
stellar bulges through two different processes: a) gas inflows 
triggered during interactions at the ODPs, and b) collision of 
baryonic cores, once the ODPs are over. The first contribute 
to the secular evolution of bulges and have been found to 
be more important at higher z, as previously discussed (see 
Figs. 6; note that these inflows are episodic). Together with 
this secular component, there is a coUisional contribution 
produced during the actual mergers of the main baryonic 
clumps, so that after a merger event, the stellar population 
of the GLO involved in the merger consists of the contribu- 
tion of both components, plus the old stars. The relative im- 
portance of each of these stellar components would depend 
on the previous evolutionary paths of the objects and their 
physical internal characteristics at the time of each merger, 
particularly on those of their potential wells (see subsection 
5.1). 

According to this scenario, the stellar bulges will have 
stellar populations with different ages. For the particular 
models analyzed in this paper, 50% of the stars in the bulge 
arc formed at z > 0.5. If we assume that, morphologically, 
the bulge component is at place when 50% of its total stellar 
mass is formed, then z « 0.5 will be its average age in these 
CDM scenarios. However, if the older stellar populations 
are taken to be age indicators, then the time of formation is 
displaced to much higher z. 

The time of formation of the bulge component could be 
affected by the cosmological model as well as by the treat- 
ment of supernova effects. In particular, the lower normal- 
ization parameters (ag) adopted for these simulations con- 
tribute to produce a late formation of bulges. In a higher 
normalization scenario, mergers would occur at higher z, 
with an earlier formation of central mass concentrations. 
Note also that higher numerical resolution can push the age 



towards fiighcr rodshift. Hence, our age indicators should be 
taken as lower limits. Nevertheless, the physical processes 
and, principally, the fact that stellar bulges are found to be 
formed naturally as the product of collapse, collisions and 
secular evolution, remain valid. Although previous works 
proposed such mechanisms to work together in the manu- 
facturing of stellar bulges, the results of the simulations re- 
ported here constitute the first proof that these mechanisms 
are actually at work in a cosmological framework and, also, 
they give us a preliminary information on how they do take 
place. The study of bulge formation by means of numerical 
simulations is very convenient, since it can provide hints on 
the physical processes involved and their interplay, while in 
semi-analytical models they have to be modelled following 
recipes known in advance (e.g., van den Bosch 2001). 

Observationally, it has been recently shown by Ellis 
et al. (2001) that at z < 0.6 the bulges of spiral galaxies 
arc bluer than ellipticals, suggesting that the stellar pop- 
ulation has been rejuvenated, specially at intermediate z. 
This has been suggested to be in contradiction with pre- 
dictions of some hierarchical clustering scenarios made by 
semi-analytical models, where spiral bulges and ellipticals 
are found to have very similar characteristics, as, in these 
models, their stellar populations have formed mostly at col- 
lapse time. However, our results show that hierarchical clus- 
tering scenarios can be reconciled with these observations, 
since the bulge stellar populations formed through mergers 
events, might be the rejuvenating agent after core collapse. 



6 DISCUSSION AND CONCLUSIONS 

The main advantage of the method reported in this paper to 
study mergers, over pre-prepared simulations, is that we an- 
alyze mergers that occur as the result of the self-consistent 
formation of the structure within a cosmological model. 
Hence, the physical properties of the colliding objects and 
the parameters of the encounters are those that correspond 
to galactic objects evolving in the adopted framework, with- 
out any external manipulation. 

In particular, the processes of gas infall at scales of hun- 
dreds of kpcs and g£is inflows in discs, at smaller scales, ap- 
pear in a natural way in our simulations. Both of these pro- 
cesses axe a consequence of evolution at larger scales, from 
initial conditions that sample a given spectrum of primordial 
fluctuations (CDM in this case) at very high z. The impor- 
tance of these processes in the SFR history of GLOs lies in 
the role they play, not only as suppliers of gas to be turned 
into stars in GLOs, but also in the fact that they provide 
the baryon mass to form a central compact concentration, 
that, later on, makes it difficult any further gas inflow. It 
is this possibility of analyzing different phenomena involved 
in galaxy formation, and the dynamical and hydrodynam- 
ical connections among them, that renders attractive the 
self-consistent approach to study mergers reported in this 
article. 

Our results suggest that the internal structure of the 
GLOs affects its SFR history, as reported by MH96 among 
others. Because a disc structure is a natural attractor for 
the dissipative component, its stability is crucial for its later 
evolution. The possible formation of radial inflows which can 
re-distribute the gas within the central regions and accumu- 
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late it in very short time-scales, directly aflFects the SFR in 
the systems. 

In agreement with previous works, wc find that some 
mergers can induce star formation in two stages, depending 
on the internal mass distributions of the parent galaxy. The 
second burst is simultaneous with the collision of the bary- 
onic cores, while the first one occurs during the orbital decay 
phase of the satellite. Wo also found that the starburst effi- 
ciency of the primary bursts tends to be larger than that of 
the secondary ones, implying that induced starbursts during 
the orbital decay phase can be more efficient at transforming 
gas into stars than those produced in the actual collision of 
the baryonic cores. Therefore, although the rate of gas cool- 
ing and accretion onto the main baryonic clump, at scales of 
say ~ 100 kpc, can be mainly determined by the dark matter 
halo (White & Frcnk 1991), how this gas distributes itself 
and responses to external and/or internal factors within the 
central GLO regions is a small scale process and it is affected 
by its particular evolutionary history (DTTS98; SalOl). 

A specific result of our method of analysis is the finding 
that a self-regulated process among gas, stars and dark mat- 
ter components is settled. This self-regulation is a basic clue 
to understand bulge formation that can only be studied in 
a self-consistent scenario. Our results show that single SBs 
tend to occur in systems with larger M^tai/M^t^^, which is 
pointing out the fact that these systems have suffered gas 
inflows in previous phases of evolution which have fed their 
bulges. As a consequence, the dark matter has also got more 
concentrated as the response to the baryon concentration 
(see Tissera & Domingucz-Tcnrciro 1998). Conversely, sys- 
tems with shallower potential wells tend to be in the pro- 
cess of forming the bulge. During this process, the whole 
potential well gets deeper. The three components, gas, stars 
and DM, work together. It is the response to their coupled 
evolution that can be correlated with the effects a merger 
event may cause on the dynamics of the systems and, as a 
consequence, on their stax formation histories. Other char- 
acteristics of mergers, such as orbital orientation, may also 
have influence on the response of the systems as reported by 
Barnes & Hernquist (1996). Their detailed study within a 
hierarchical scenario will be carried out in a separate paper. 

Wc found that in a hierarchical clustering scenario, 
there is a natural evolutionary sequence: mergers occurring 
at higher z when the systems are not that concentrated and 
stellar bulges not well-formed, lead to gas inflows that trig- 
ger SBs and feed the first bulges. As a consequence, the total 
potential wells concentrate. When the objects reach lower 
z, they tend to be more stable because of their important 
stellar bulges and concentrated potential wells that prevent 
further gas inflows. At this stage, mergers also tend to be 
minor. This evolutionary sequence suggests that starbursts 
induced during the ODP are likely to be more common at 
higher redshifts: and because of their higher efficiency in 
transforming gas into stars, we predict that their impact on 
galaxy formation could be very important, particularly, on 
the formation of spiral bulges. 

Finally, our results suggest that hierarchical clustering 
scenarios can be reconciled with recent observational results 
of high redshift spiral bulges and ellipticals in the field, since, 
in such scenarios, stellar bulges would form by the contri- 
bution of collapse, collisions and secular evolution, each one 



with a relative importance determined by the evolutionary 
path of each system. 
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Table 1. Numerical Simulations 



s 


CI 


0-8 


c 


S.l 


A 


0.40 


0.01 


S.2 


A 


0.40 


0.10 


S.3 


B 


0.67 


0.01 


S.4 


C 


0.67 


0.01 



Table 2. Main Characteristics of Starbursts 



GLO 


JVdark 


JVbar 


Cstar 


■'Wburst 


"njurst 


Msat/Mpro 


Mstar/Mbar 


A^AarZ-'^star" 


^*burst 


^burst 


1 


9181 


1335 


18.75 


0.55 


4.36 


0.16 


0.12 


0.20 




0.50 








33.01 


4.06 


19.53 


0.23 


0.10 


0.30 


27.09 


0.41 








15.47 


1.12 


8.61 


D 








0.15 


2 


7310 


1059 


24.76 


0.96 


3.53 


1.08 


0.14 


0.10 


5.42 


0.63 








16.52 


0.10 


1.99 


D 








0.54 








46.43 


4.63 


17.60 


0.39 


0.26 


0.35 




0.47 








12.38 


0.26 


5.76 


D 








0.29 








17.53 


1.25 


10.11 


0.11 


0.28 


0.90 




0.08 


3* 


6317 


917 


13.41 


0.91 


11.50 


0.40 


0.22 


0.40 


16.67 


0.56 








28.89 


2.18 


11.21 


D 








0.34 


4* 


7115 


908 


29.92 


1.44 


13.50 


0.16 


0.24 


0.60 




0.34 








12.38 


0.65 


10.39 


0.16 


0.27 


0.85 




0.14 


5* 


5215 


739 


8.25 


0.49 


2.48 


0.43 


0.27 


0.42 


25.01 


0.41 








19.60 


0.89 


12.4 


D 








0.34 








22.70 


0.91 


5.14 


0.22 


0.21 


0.75 




0.19 


6* 


6184 


798 


16.50 


1.71 


11.50 


0.40 


0.15 


0.25 


2.58 


0.78 








17.53 


1.09 


8.11 


D 








0.52 


1 


9199 


1307 


42.30 


0.68 


3.94 


0.18 


0.58 


0.20 




0.56 








43.33 


2.76 


13.40 


0.22 


0.61 


0.63 




0.43 


2 


7383 


1071 


52.62 


1.56 


8.86 


1.06 


0.51 


0.15 


13.01 


0.81 








30.96 


1.72 


6.01 


D 








0.56 








71.19 


4.18 


19.08 


0.38 


0.54 


0.40 




0.39 








66.93 


0.36 


8.13 


0.11 


0.59 


0.95 




0.09 


3 


6501 


918 


30.95 


1.56 


11.21 


0.43 


0.62 


0.90 




0.34 


4 


7144 


921 


39.21 


1.87 


7.47 


0.10 


0.66 


0.55 




0.37 








41.12 


1.53 


13.99 


0.21 


0.46 


0.88 




0.13 


5 


5294 


744 


40.24 


2.11 


15.03 


0.51 


0.78 


0.45 


15.80 


0.49 








33.02 


1.43 


19.30 


D 








0.30 


1 


3688 


550 


23.73 


0.57 


0.81 


0.31 


0.10 


0.30 


37.00 


1.47 








5.16 


0.05 


1.23 


D 








0.58 


2 


5665 


1056 


8.25 


0.60 


14.25 


0.58 


0.15 


0.34 


31.11 


0.65 








5.16 


0.44 
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Table 3. Mean Starburst Efficiencies 

S S.l S.2 S.3 S.4 

'^stB.r/^stt? ^ 0-40 0-39 0.56 0.16 0.11 

M^^^JM^=° > 0.40 0.08 0.39 0.03 0.03 
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Figure 1. Star formation rate for a typical galaxy-like object 
experiencing a double starburst (a) and a single one (b) in simu- 
lation S.l as a function of lookback time. The spatial separation 
between the mass centres of the progenitor and satellite clumps 
are plotted {dashed lines). The solid lines represent the quiescent 
star formation component. 
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Figure 2. Number of baryonic particles witliin radial distance 
r for GLO 2 (a) in S.4 and GLO 5 (b) in S.l between Zo and zi. 

Specific angular momentum component j^^i along the direction of 
the total angular momentum of the progenitor for each baryonic 
particle within Ri 100 kpc of GLO 2 (c,e,g) and GLO 5 (d,f,h) 
versus their positions from Zo to zi. Circles: gas particles, stars: 
star particles, asterisks: counterrotating stars. 



© 0000 RAS, MNRAS 000, 000-000 



14 Tissera et al. 



CO 



b 



•if 



o 



o 



O 



O Oo 







•t; O 



OA® ▲ 







0.2 



0.4 



0.6 



0.8 



star/ 



star 



Figure 3. Ratio between the strengths of the double stcirbursts 
("'star /""star ) ^^"^ fraction of stars already formed in the pro- 
genitor (A^star/^star) ^ ^'^'^ merger. Single bursts have 
been assigned ''■star/''^star ~ ^- GLO 1 in S.3 has not been in- 
cluded for the sake of clarity I f^'^i^^^ = 4.6). GLOs in S.l: 
open circles, S.2: filled triangles S.3; filled squares and S.4; open 
stars. 
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Figure 4. Ratio between tlie strengths of the double starbursts 
("^star / "^star ) versus the circular velocity at 2 kpc of the a) total 
mass distribution (wtot), b) dark matter (i^dm) S'ld c) the gas 
('L'gas) for all bursts in S.l {open circles) and S.4 (open stars). 
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Figure 5. The fraction of gas consumed into stars during star- 
bursts, Mburst /A'/gas ■ as a function of the fraction of stars already 
in place in the progenitor at the z of the merger (^M^sftarZ-'^star)- 
Single bursts and primary components of double ones have been 
included. We find that, in a given simulation, higher burst ef- 
ficiencies tend to occur in GLOs that lack an important stellar 
mass concentration (see Fig. 3 for feature code). 
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Figure 6. a) Msat/Mpro and b) M^^^JM^^;^^ as a function of the 
•^burst fof single and primary components of double ones (symbols 
with arrows). See Fig. 3 for feature description. 
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